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Abstract 

We study the latent heat of the liquid-gas phase transition in symmetric nuclear matter using 
self-consistent mean-field calculations with a few Skyrme forces. The temperature dependence of 
the latent heat is rather independent of the mean-field parametrization and it can be characterized 
by a few parameters. At low temperatures, the latent heat tends to the saturation energy. Near 
the critical point, the latent heat goes to zero with a well-determined mean-field critical exponent. 
A maximum value of the latent heat in the range I ~ 25 — 30 MeV is found at intermediate 
temperatures, which might have experimental relevance. All these features can be explained from 
very basic principles. 

PACS numbers: 21.60.Jz, 21.65.-f, 64.60.F- 
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I. INTRODUCTION 



The study of the liquid-gas phase transition in homogeneous symmetric nuclear matter 
provides interesting links between statistical mechanics, quantum physics and nuclear dy- 
namics [l;3]. The liquid-gas phase transition picture can, in principle, be applied to nuclear 
systems 4|. The nucleon-nucleon (NN) interaction has a qualitative similar structure to the 
inter-atomic one: it is repulsive at short relative distances and attractive at intermediate 
and long distances [5|. At zero temperature, this structure drives the particles to stay at 
a given distance from each other, thus leading to a structured liquid phase. As tempera- 
ture increases, thermal fluctuations overcome the interaction effects and the system tends 
to gasify, thus causing a phase transition jfj]. 

The only way to study thermodynamical properties of nuclear systems on Earth is via 
intermediate and heavy ion collisions jl[ Q], which unfortunately are particularly fast (~ 
10 -21 s) processes and very difficult to model. Unlike in other branches of physics, one cannot 
define or prepare samples for a thermodynamical experiment. Furthermore, in a collision 
the state of the system changes continuously in time, and it is difficult to establish whether 
the system has actually reached a state of equilibrium or not 7|, |8| . Finally, the fact that we 
are dealing with finite nuclear matter causes additional difficulties on the thermodynamical 



interpretation of the results 



9]. 



Experimentally, the challenge lies in being able to control the equilibrium, i.e. in deter- 
mining the state variables of the system, like the temperature, the pressure or the density of 



the system 
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llj . Moreover, one also needs to find suitable observables that help to iden- 



tify the liquid-gas phase transition. Candidates include, among others, the critical behavior 
of fragment partitions energy fluctuations on an event by event basi s | l5l| . nuclear 



calorimetry [16|, [17J and bimodality of the largest fragment distribution 18|, |l9j. The evi- 
dence gathered with these different experimental techniques points towards the existence of 
a liquid-gas phase transition for nuclear systems at densities below the empirical saturation 
density, po — 0.16fm -3 , and temperatures around 5 — 10 MeV. 

From the theoretical side, symmetric nuclear matter at finite temperature provides a 
first qualitative picture of the thermodynamics of nuclear systems and, particularly, of the 
liquid-gas phase transition 

a a. 

Because of its relative simplicity, nuclear matter has 
been the subject of numerous investigations. From a microscopic perspective, different 
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many-body techniques (including Dirac-Brueckner 20] and Brueckner-Hartree-Fock 2 If . 



Self- Consistent Green's functions 122 



23] or variational calculations 24|) have been used 



to compute the thermodynamical properties of nuclear matter starting from basic nucleonic 
degrees of freedom and phase-shift equivalent realistic NN interactions. So far, very different 
predictions for critical properties have been obtained using different approaches or even 
different interactions within the same approach 22J . A more efficient way to get insight in the 
thermal behaviour of nuclear matter is by means of effective interactions in the framework of 
a Hartree-Fock approximation, either relativistic or non-relativistic (26 , [27 ] . Mean-field 
calculations are much faster to implement numerically and give access to a wider range of 
phenomena (critical exponents, for instance). Moreover, the mean- field parametrizations are 



directly connected to nuclear structure 



28| . Recently, one of us has studied systematically 



the dependence on the effective interaction of the liquid-gas phase transition of nuclear 
matter with Skyrme and Gogny mean-fields 291 ]. 

In the context of homogeneous nuclear matter, one generally assumes that thermody- 
namical equilibrium takes place between two infinite pieces of nuclear matter, one of them 
belonging to the liquid phase and the other to a gaseous state. This is quite a crude approxi- 
mation that can only provide an average, qualitative understanding of the thermal properties 
of the liquid-gas phase transition. As a matter of fact, one should consider that the low- 



density phase of dense nuclear matter is not necessarily homogeneous [30J. Its description 
should allow also for the formation of light clusters of nucleons, including deuterons, tri- 
tons, helions ( 3 He) or a-particles 3lH34|. Moreover, finite size and Coulomb effects are very 
important to determine the thermal properties and the critical behavior of finite fermionic 
systems [35]. 



Very little attention has been paid so far to the study of the latent heat of nuclear matter, 
which is the main goal of this work. This is a very interesting quantity, since it can provide a 
further characterization of the liquid-gas phase transition, and it can be potentially extracted 
both from experiments and theory. Experimentally, the latent heat might be, in principle, 
read out from the length of the plateau in the caloric curve (i.e. in the temperature vs. 
excitation energy curve). This would suggest values of I ~ 4 — 8 MeV 



16 



17, 



36|. Note, 



however, that this p 
transition picture 



ateau can also be explained in other terms, rather than in a phase 



37H39|. Alternatively, the latent heat can also be obtained from the heat 
capacity versus excitation energy, as shown in Ref. [IS] , in which case / ~ 2 — 4 MeV. 
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Recently, bimodality in the charge distribution of the heaviest fragments has been used to 



extract a latent heat value of I ~ 8 MeV 



19]. In any case, one has to be aware of the fact 



that, for the latent heat to be defined, the reaction employed to extract it should not be an 



isentropic process 



40]. 



In the case of infinite nuclear matter, the latent heat per particle accounts for the amount 
of heat needed to take a nucleon from the liquid to the gas phase. In this work, we want to 
study this quantity and provide estimates for its typical values in nuclear matter. We will, 
in particular, describe its basic temperature dependence with different Skyrme mean-field 
parametrizations. We also want to link the behavior of this quantity to the properties of the 
underlying NN interaction. Ultimately, our aim is to understand the behavior of the latent 
heat from basic quantum statistical mechanics arguments. 

In the next section, we discuss very briefly the formalism needed for the calculation of 
the latent heat within the Skyrme-Hartree-Fock approximation. After that, in Section III, 
we present and discuss the self-consistent results, which is followed by a section devoted to 
the analysis of the low temperature and critical behavior of the latent heat. A summary 
and final conclusions are presented in the last section. 

II. FORMALISM 

To describe efficiently finite nuclei and nuclear matter at the same time, one generally 
relies on phenomenological interactions, adequate for the Hartee-Fock (HF) approximation. 
The effective 
in the 1950's 



orces used in the following are of the Skyrme type. They were introduced 



41 



42J, and have been intensively used in the literature |43|]. Alternatively, 



one can obtain equivalent results by formulating the problem in terms of density functionals 
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44]. 



In this work we deal with a Skyrme effective interaction: 

+ t 2 (l + x 2 P„)t '6(l>)t + Wl + x 3 P a )[p(R)] a S(-p'), (1) 

o 

where ~Pt = (rt + Tj)/2, ~P* = rt — rj, k = (Vj — Vj)/2i is acting on the right; 
k = — — ^~j)/2i, acting on the left; P a = (1 + a% ■ a%)/2 is the spin-exchange op- 
erator (<r are Pauli matrices;), and p = N/Q, the total baryonic density. The parameters 



to,ti,t 2 ,ts,xo,xi,X2,xs,a are numerical constants fitted to reproduce, in general, the sat- 
uration properties of nuclear matter and structure properties of closed shell nuclei j^ . 
Three-body interactions are effectively accounted for by the last density-dependent term. 
Mean-field calculations with effective Skyrme forces reproduce in a satisfactory way the 
structure of a wide range of nuclei [281 ] . 

Due to the translational invariance of uniform nuclear matter, single-particle (s.p.) states 
are appropriately described by plane waves, and the relevant quantum numbers are the s.p. 
momentum, k, as well as the spin and isospin projections. At finite temperature within the 
Hartree-Fock approximation, momentum states are occupied according to the Fermi-Dirac 
distribution, 

n(k,T) = \ , (2) 

where T is the temperature of the system; /i, its chemical potential, and e(k), the s.p. en- 
ergy. Including the rearrangement term arising from the density dependence of the effective 
interaction, the latter is written as: 



t + ^(a + 2)t 3 p c 



(3) 



where the effective mass, 



m 



m 1 + §^ (3ti + 5t 2 + 4t 2 x 2 ) p 
and the kinetic energy density, 



(4) 



T= j£r J d3kk2n ^ T) > (5) 

have been introduced. The integral of n(k, T) over the available phase space at finite tem- 
perature T gives the total density, p: 

P = j£y J d 3 kn(k,T), (6) 

where v is the degeneracy [y — 4 for spin and isospin saturated nuclear matter). This 
condition determines the chemical potential at a fixed external density. The calculation of 
the single-particle energy and the density normalization condition defines a self-consistent 
process. Note, however, that due to the simple structure of the Skyrme interaction, self- 
consistency at the HF level is already achieved at the first iteration. 



Having fixed the chemical potential, the momentum distribution and the kinetic energy 
density are fixed. From these quantities, one can immediately compute the total HF energy 
per particle: 



e( P ,T) 



h 2 T 3 

2m* p 8 



P 
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hp 



(7) 



and the entropy per particle: 



d 3 k 



n(k, T) \nn{k, T) + (1 - n(k, T)) ln(l - n(k, T)) 



(8) 



(2?r) 3 p 

In turn, these give access to the free energy, / = e — Ts, from which all the remaining 
thermodynamical properties can be computed. In addition to the chemical potential, /i, 
the pressure, P, is also necessary to study the liquid-gas phase transition. One should take 
into account that this procedure is thermodynamically consistent, in the sense that the 
chemical potential extracted from the normalization of the density, Eq. OH]), coincides with 
the chemical potential derived from the free energy, p = f + p (df /dp), provided that the 
rearrangement term is properly included. 
The Clausius-Clapeyron formula, 



g Pi) \ dT J, 



(9) 



gives access to the latent heat per particle 



45]. This is expressed in terms of the product 



of the temperature, the difference between volume per particle of the two phases and the 
derivative of the pressure with respect to the temperature along the coexistence curve. 
Alternatively, the latent heat can also be computed as the amount of heat that is needed, 
at a fixed temperature, to transfer one nucleon from the ordered (liquid) to the disordered 
(gas) phase. Since this process happens at constant chemical potential and pressure, the 
heat change only involves the difference in entropy per particle of the two phases: 



I = T(s Q - s{] 



(10) 



For nuclear matter calculations, this formula is numerically more stable than the Clausius- 
Clapeyron one. Values of the vaporization specific latent heat for common liquids and gases 
are in the region of 10 — 5000 kJ/kg when measured at the normal boiling point. In contrast, 
as we shall see in the following, nuclear matter has a maximum specific latent heat of the 
order of ~ 30 MeV, i.e. ~ 3 x 10 12 kJ/kg, which is orders of magnitude higher and among the 



highest in nature (an exception being the quark-gluon plasma). The origin of this extremely 
large value can be traced back to the strong force, which binds the nucleons tightly. 

A qualitative model that justifies these relatively high values and provides insight on the 
temperature dependence of the latent heat can be obtained using basic ingredients. From 
Eq. (HO I) , we see that the latent heat can be computed from the difference of entropies times 
the temperature. On the one hand, let us take a classical approximation for the entropy of 
the gas, 

This approximation should be valid because the gas coexistence density is, in general, quite 
low, which is also why interaction effects are neglected. On the other hand, the liquid phase 
is closer to saturation density and one could use the Sommerfeld expansion to compute its 
thermodynamical properties. The entropy in the liquid phase would then yield: 

7T 2 T 

s i = TT> (12) 

h 2 k 2 

where Eq = -^f-. We have ignored interaction effects (which would appear in the form of 
an effective mass) and we take kp at a reference density equal to saturation density. Within 
this approximation, the latent heat becomes: 

/ = - M (T) + \T - |V . (13) 

The chemical potential is the same for the liquid and the gas phases, due to the coexistence 
conditions (see Eq. (fl8l) below). One can therefore use the Sommerfeld expansion in the 
liquid branch to find the corresponding chemical potential: 



A*z — £o 



^2 / rp \ 2' 
71 I 1 x 



1 12 U 



U . (14) 



Note that, again, we fix all densities to be at saturation and that Uq is the attractive 
contribution needed so that /^(T = 0) = e(p ,T = 0) = — e . Using Eq. f[T4|) into Eq. ffTB"]) , 
one finds the following expression for the latent heat: 

'- + ^-^' (15) 

Within this very crude model, the latent heat is a quadratic function of the temperature. 
Loosely speaking, classical effects (i.e. those due to Eq. (Till ) tend to increase the latent heat 



linearly in temperature. Thermal effects arising from the degenerate expansion compensate 
this term and lead to a maximum, occurring at a temperature, 

T max = ^ ~HMeV. (16) 

The maximum value of the latent heat is then given by the sum of two terms: 

/ max = e + -^r ~30MeV. (17) 

In addition to the term due to the saturation energy, an equally important term (in sign and 
size) appears as a result of the competition between classical and degeneracy effects. This 
increases the value of the maximum latent heat substantially, up to around 30 MeV. In the 
following, we will see that the predictions of this simple model are well fulfilled by Skyrme 
mean-field calculations of the liquid-gas phase transition. 



III. RESULTS 



The liquid-gas phase coexistence and the critical properties of nuclear matter studied in 
a mean-field approximation with Skyrme forces are well-known 29] . Since the emphasis 
here is on the latent heat itself, we will only consider a few representative effective forces. 
We shall also briefly comment on some well known results on the phase coexistence before 
considering the latent heat. These will be helpful for our analysis. As a first example, we 
choose the Skyrme force BSkl7 46[, which gives a very accurate description of the masses 
of nuclei all across the mass table. Four pressure isotherms calculated with BSkl7 are 
reported in the left panel of FigJTJ These include the representative T = (solid line) and 
the critical T = T c (dashed line) isotherms. Two more isotherms, one below (dotted line) 
and one above (dash-dotted line) the critical one, are also displayed. The isotherms show 
the well known shape associated to a liquid-gas phase transition. At T = 0, the liquid at 
saturation density has zero pressure and, therfore, it is in equilibrium with a zero density 
gas. As the temperature rises, the gas coexistence density shifts to finite values and the 
coexistence region shrinks until the critical temperature is reached. Below T c , all isotherms 
present a mechanically unstable region, where the pressure decreases with density. At T c , 
phase coexistence is not possible anymore and the system vaporizes completely. In terms 
of isotherms, one finds that, for T > T c , the pressure becomes a monotonically increasing 
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BSkl7 Isotherms Reduced Isotherms 




0.05 0.1 0.15 0.5 1 1.5 2 2.5 



p (fm 3 ) P ; Pc 

FIG. 1: (Color online) Left panel: four different isotherms for the BSkl7 Skyrme force, T c = 
15.6MeV. Right panel: reduced isotherms at T = T c and T = 0.7 T c for the BSkl7 (solid lines) 
and ZR1 (dash-dotted lines) forces. Isotherms of the Van der Waals model (dotted lines) are also 
shown for comparison. 



function of density. For BSkl7, the critical temperature turns out to be 15.6 MeV, similar 



to the T c of a wide range of modern Skyrme forces 



29]. 



The behaviour of these isotherms is reminiscent of the Van der Waals equation of state 
(EoS) for real gases which, unlike the ideal gas EoS, takes into account the non-zero size of 
the molecules of the gas [6|]. The excluded volume has a repulsive effect, in contrast to the 
attractive inter-molecular force. This causes the formation of regions of instability which 
result in a liquid-gas phase transition. To appreciate the similitudes and differences between 
the Van der Waals and the self-consistent mean- field EoS, we plot in the right panel of Fig. [1] 
the reduced isotherms (P/P c vs. p/p c ) at T = 0.7 T c and at T c . In addition to the Van der 



Waals case, we consider the EoS obtained with BSkl7 and with the ZR1 



471 ] forces. The 



latter has the highest critical temperature (T c = 22.98 MeV) among the large set of effective 
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BSkl7: liquid-gas phase coexistence 




P (fm ") 



FIG. 2: (Color online) Coexistence line (solid) and spinodal line (dotted) obtained from the BSkl7 
force. 



forces analyzed in Ref. 29]. 

By studying the pressure in reduced units (P/P c , instead of P) we expect to highlight 
possible resemblances between several EoS. The Van der Waals model describes the EoS of 
monatomic gases and liquids within a reasonable distance above and below their critical point 
js]. Moreover, according to the principle of corresponding states, if we measure pressure, 
volume and temperature in units of P, p c and T c , the EoS becomes universal, i.e. it is the 
same for a wide range of substances |6j. One might wonder whether there is something like 
a principle of corresponding states for EoS derived from different mean-fields. The right 
panel of Fig. [TJ provides an insight into this matter. For T = 0.7 T c , beyond the gas phase 
there is a relatively important disagreement between all the reduced EoS. It is only close to 
the critical point that the EoS coincide. This is a consequence, as we will see later, of the 
fact that the Van der Waals and the self-consistent mean-field models have the same critical 
exponents [29 ]. 



When nuclear matter is heated up, a phase coexistence develops between a relatively 
high-density liquid phase and a low-density gas. In order for such equilibrium to exist, the 
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chemical potential and the pressure of the two phases should be equal: 



p(p g ,T) = »( Pl ,T) , P( Pg ,T) = P{ Ph T) . (18) 

At a given temperature, T, the solution of this set of equations provides a density couple, 
(p g ,pi), which defines the coexistence densities of the two phases. In Fig. [21 we plot the 
coexistence phase diagram (solid line) in the (p, T) plane for the BSkl7 force. As it has 
already been mentioned, at T = the liquid at saturation density is in coexistence with 
a zero-density gas. The gas coexistence density grows and the liquid coexistence density 
decreases as temperature rises. The densities of the gas and liquid phases join at the critical 
point, above which the liquid-gas phase coexistence disappears. Together with the coexis- 
tence line, we also plot the spinodal line (dotted line), which marks the boundary between 
thermodynamically stable and unstable states of matter. 

The latent heat corresponding to the liquid-gas phase transition obtained with BSkl7 
using Eq. ( flOl) is reported in the left panel of Fig. [3] (solid line). One can observe that the 
latent heat has a characteristic bump shape as a function of temperature. In the T = limit, 
/ has a finite value. For BSkl7 (solid line), it grows with temperature up to a well defined 
maximum at T = 8.7 MeV (with Z max = 29.9 MeV) and then sharply goes to zero at T = T c . 
To appreciate better the dependence of the latent heat of nuclear matter on the different 
Skyrme parameterizations, we have also performed calculations of the coexistence line and 



latent heat for various other forces. In addition to t 
we have chosen two more interactions: SLy9 



48 



re already mentioned BSkl7 and ZR1, 



491 ]. which incorporates, by construction, 
the behavior of a microscopically derived EoS of neutron matter [|3(3] and has a very low 



critical temperature 29j; and LNS, that reproduces the properties of BHF calculations of 
nuclear matter [511 ]. The results obtained for their latent heats are also reported in the left 
panel of Fig. [3] Note that all forces produce a latent heat with a similar qualitative behavior 
that, as we shall discuss below, can be understood from basic principles. 

Let us start the discussion with the zero temperature limit of the latent heat. As men- 
tioned previously, / measures the heat that needs to be provided to the fluid to transfer a 
nucleon from the liquid to the gas phase. Since the equilibrium coexistence density of the 
gas asymptotically goes to zero when taking the T — Y limit, I becomes the amount of 
heat needed to extract a particle from the system in that limit. In other words, for T — > 0, 
/ —7- —jjL. At the saturation density of nuclear matter, p , one actually finds p = — e and, 

11 



Latent heat Reduced latent heat 




FIG. 3: (Color online) Left panel: latent heat using four different Skyrme force parametrizations. 
Right panel: latent heat in reduced dimensionless units using the same forces. 



consequently, the latent heat is equal in absolute value to the saturation energy, I = eo, 
when getting close to T = 0. This result is independent of the interaction or the many-body 
approximation used to describe nuclear matter, as long as thermodynamical consistency is 
fulfilled. This suggests that, if one wants to obtain interaction- independent results for the 
latent heat, I might be normalized to eo (see right panel of Fig. [3]). 

We observe that for all Skyrme parametrizations the qualitative behaviour of the latent 
heat is very similar (left panel of Fig. [3]). The latent heat matches the value of the binding 
energy at T = 0, then rises for small temperatures. The initial rise is linear and the slope 
seems to be independent of the Skyrme parametrization. As we shall see below, the slope of 
I close to T = is a model-independent feature that can be understood from fundamental 
arguments. Further up in temperature, I reaches a maximum and then drops to zero at 
the critical point, where the difference between the liquid and gas phases disappears. The 
results with this limited set of forces seem to indicate that the position and magnitude 
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of the maximum in I depend on the specific value of the critical temperature. Broadly 
speaking, higher values of T c shift the position and height of the maximum to larger values. 
This is particularly clear for the case of the ZR1 mean-field, which has the largest critical 
temperature. A calculation with a wider set of mean-fields, not shown here for simplicity, 
confirms this tendency. 

A plot of the latent heat in reduced units, l/eo, as a function of the reduced temperature, 
T/T c , is presented in the right panel of Fig. |3j Similar to the right panel of Fig. [H the 
reduced plot is helpful in highlighting the dependence of the latent heat on the different 
EoS. In general terms, we observe that the large dependence on the mean-field is eliminated 
to a large extent in the dimensionless plot. Close to T = 0, the linear slope of I is changed 
due to the fact that different mean- fields have different saturation energies, e$. Yet, near 
the maximum, the reduced latent heats show a much smaller deviation compared to the 
absolute ones. For all Skyrme forces, the latent heat tends to peak within a limited region 
of temperatures, T/T c ~ 0.5 — 0.6. Moreover, the peak is also quite narrowly distributed 
around the value l m ax/^o ~ 1.7 — 2, which suggests that the latent heat is more determined 
by thermal correlations than by effective forces. Finally, as the temperature reaches the 
critical value, the latent heat falls to zero with a very similar temperature dependence for all 
forces. This identical behavior can be explained in terms of critical exponents [see Eq. (13ip ]. 
As it will be shown later, within the mean-field approximation, the critical exponents of all 
latent heats are the same close to the critical point. 

As mentioned in the previous section, the Clausius-Clapeyron formula requires the evalu- 
ation of the derivative of the pressure with respect to the temperature along the coexistence 
line. The pressure along coexistence is commonly referred to as the vapor pressure. Its be- 
havior in reduced dimensionless units, P v /P c , along the phase transition is shown in Fig. H] 
for the four mentioned Skyrme parametrizations. The vapor pressure is a well behaved func- 
tion of the temperature that grows from zero to P c as temperature increases from zero to T c . 
For low temperatures, P v rises very slowly and at T = 0.5 T c it is only 10 % of P c . Above 
this temperature, a steady increase brings the vapor pressure very rapidly up to P c . Note 
that the last portion of this increase is basically linear, in accordance to Eq. ( 130]) below. 
We have checked the numerical and thermodynamical consistency between the values of the 
latent heat obtained with the Clausius-Clapeyron equation (Eq. (jUJ)) and those given by the 
difference of entropies of the gas and the liquid phase (Eq. (fTOl) ). 
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FIG. 4: (Color online) Vapor pressure versus temperature in reduced dimensionless units obtained 
from the equilibrium conditions Eq. (|18p . 



A basic ingredient in the evaluation of the latent heat is the difference between the 
entropy of the gas and that of the liquid phase [see Eq. (flQj) ]. This difference, which has 
to be evaluated along the coexistence line, is shown in Fig. [5] for the BSkl7 interaction. 
Entropies are plotted in a logarithmic scale. The difference between the entropy of the gas 
and that of the liquid is always a positive quantity that decreases with temperature and goes 
to zero as the system reaches the critical temperature. The low-temperature limit of this 
difference can be studied analytically and will be discussed in the next section. Note that, 
within this limit, the difference between the gas and liquid entropy is largely dominated by 
the entropy of the gas. Also notice that the variation of the liquid entropy, which should 
be zero at T = 0, is much smaller than the variation of the entropy associated to the gas 
phase (note the logarithmic scale). While the liquid entropy along the coexistence line is 
an increasing function of temperature, the gas entropy is a decreasing one. Their difference, 
however, is dominated by the gas entropy and becomes a decreasing function of temperature. 
As a matter of fact, in the context of classical gases, it is customary to neglect the liquid 
contribution, since its entropy (or, in terms of Eq. fl9]), its inverse volume) is negligible with 
respect to the gas one |45j. In nuclear matter, this approximation would only be valid up to 
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Entropy per nucleon 




T (MeV) 

FIG. 5: (Color online) Logarithmic plot of the difference of gas and liquid entropies (solid line), 
entropy of the liquid (dashed line) and entropy of the gas (dotted line) along the coexistence line 
using the BSkl7 Skyrme force. 

T/T c ~ 1/3, which would lead to a maximum error in the calculation of the latent heat of 
about 10%. For higher temperatures the error induced by this approximation would already 
reach the 30% for the maximum of the latent heat. Finally, let us stress that the maximum 
that appears in the latent heat at intermediate temperatures is a subtle result, arising from 
the product of the (linearly increasing) temperature times the (decreasing function of T) 
difference of entropies. 

IV. ANALYTICAL LIMITS OF THE LATENT HEAT 
A. Low temperature behaviour 

As we have already discussed, all Skyrme forces produce a qualitative similar behaviour 
of the latent heat. The value in the zero temperature limit (/ = eo) is well understood from 
basic arguments. In this limit, the liquid phase is at the saturation density and therefore has 
zero pressure. The heat needed to transfer a nucleon from the homogeneous liquid nuclear 
matter phase to the zero density gas (vacuum) is just the chemical potential of the liquid 
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which, at saturation, coincides with the binding energy per particle. In this subsection we 
will show analytically that this intuitive value is the result of a delicate balance. 

In our derivation, we will not only evaluate the value of I at T = 0, but also its derivative 
as a function of temperature at this point. This derivative is basic in understanding the 
existence of a maximum in I. The argument goes as follows: at T = 0, the latent heat is 
finite and positive. As T — > T c , however, the latent heat must go to zero, since the liquid-gas 
phase transition disappears. The existence of the maximum is therefore necessarily related 
to the way / departs from its T = value. If the slope at T = is positive, the function will 
first increase with T and, in order to reach the zero value at T c , at least one maximum will 
have to develop at some intermediate temperature. We shall also show that the slope of I 
at T = is not only positive, but model-independent. 

Let us consider the limit of low temperatures of the latent heat. As shown in Fig. [2J as 
T approaches to zero, the dense liquid phase is in equilibrium with a very low-density gas. 
Due to its diluteness, we shall assume that the gas is in the classical regime and interactions 
are not relevant anymore. The thermodynamical properties of the gas are then given by 
ideal gas expressions p and they can be computed analytically for a given gas density, p g , 
and temperature, T. The pressure is simply given by 



while the gas density can be obtained in terms of its chemical potential by using the following 
expression: 



If we substitute the previous expression in the pressure, we find the vapor pressure of nuclear 
matter: 



Notice that the information on the interactions in nuclear matter is contained in the chemical 
potential, fi g , which is the same for the liquid and the gas according to the equilibrium 
relation, Eq. (fl8l) . In the liquid branch, as T approaches zero, the chemical potential tends 
to the energy per particle, and one has fi g = \i\ — > — e . 

According to the Clausius-Clapeyron formula, Eq. to calculate / we need the tem- 
perature derivative of the vapor pressure along the coexistence curve. Using Eq. ( 12T|) . this 



P = PgT, 



(19) 




(20) 




(21) 
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derivative becomes: 



dP v 



m \ 3 /2 



dT 



Pg_ 
T 



T ,' 5 d % 
2 dT 



(22) 



All temperature derivatives are to be taken as derivatives along the coexistence line. Because 
of the equilibrium condition in Eq. f fl8|) . the chemical potential can be computed from 
the liquid one which, close to saturation, should be that of a degenerate Fermi gas. The 
Sommerfeld expansion then guarantees that the temperature dependence of pi is quadratic 
in T 52]. Consequently, the temperature derivative in the last term ^ ~ 0{T) and thus 
it can be neglected in the following considerations. Note that the low temperature limit of 
the previous expression, Eq. fT22|) . is zero. However, when the latent heat is considered in 
the zero temperature limit, we need to explicitly take into account the prefactors, whose 
cancelation leads to: 



lim l(T) = lim T I 



dT 



lim T- 



1 dP, 



-p g = e • 



(23) 



t->o p g dT 

When taking the limit, we have considered that the term containing the liquid density goes 
to zero because the gas density present in the derivative of the vapor pressure (see Eq. ( 1221) ) 
goes to zero in this limit, and that the linear terms in temperature within the brackets in 
Eq. fT22|) are subleading. 

Let us now compute the derivative of I close to zero temperature: 



dl 
dT 



1 



1 

Pi 



dP„ 



dT 



V - + T 



1 

Pg 



1 

Pi 



d 2 P v 
dT 2 



T 



dP v d 
lTdT 



1 

Pi 



(24) 



The T — > limit of the first term is obtained from Eq. (|22|) . The second term involves the 
second derivative of the vapor pressure, which is also easily computed from Eq. (T22]) and 
yields: 



1 

Pg 



1 

Pi 



d l P v 
dT 2 



1 

Ji2 



(25) 



The third term involves derivatives of the gas density, Eq. ( 1201) . while the liquid density is 
neglected: 



T- 



. dP i: d 



1 



Pg Pi 



J^2 



dT dT 

Collecting the different contributions, one gets : 

,2 



15, 



1 ~ — I - 3pT + —T 



dl 
lim — 

T^O dT 



T 



+ 



p_ 

J^2 



3p 15 

Y + T 



p_ 

J^2 



4/i 15 

Y + T 



5 
2 



(26) 



(27) 



17 



Latent heat fit at low T 
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FIG. 6: (Color online) Latent heat versus temperature for very low T. Symbols represent self- 
consistent mean-field results and lines their respective linear regression fits. 



The derivative of the latent heat with respect to the temperature in the limit T — » is 
therefore independent of the interaction. Moreover, it is positive and equal to 5/2. 

In Fig. O we show a numerical proof of Eq. (1271) . For the different Skyrme force 
parametrizations considered in Fig. [31 we focus on the low-temperature behaviour of the 
latent heat (symbols). The slopes for the linear regression fits (lines) of these numerical 
data are in very good agreement with the value 5/2 up to two significant digits. 

This result is not only valid regardless of the effective interaction, but it is also valid 
no matter which many-body approximation is considered. The only assumption that has 
been made is that the gas equilibrium density enters the classical regime as the temperature 
decreases. Consequently, one should get the same result in approaches that go beyond the 
HF approximation. It is also important to note that this result is independent of the system 
under study and therefore should be generically valid for the liquid-gas phase transition of 
any extended normal fermionic system. As discussed previously, the positiveness of this 
derivative necessarily implies that a maximum in the latent heat must develop. As a result, 
we expect a maximum in the latent heat of any normal fermionic system that presents a 
liquid-gas phase transition. 
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Critical exponent 
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FIG. 7: (Color online) Logarithmic plot of the latent heat versus the reduced temperature. The 
slope of the different curves gives the critical exponent. Symbols represent self-consistent mean- field 
results and lines their respective linear regression fits. 



B. Critical behaviour 



Critical exponents characterize the properties of phase transitions js), l^, 53[. Close to the 
critical point, all the thermodynamical properties can be described in terms of a handful of 
exponents. Their knowledge facilitates the understanding of the properties of systems close 
to criticality. The latent heat is not an exception and it can also be described, close to the 



critical point, in terms of a critical exponent. As shown in Ref. 45], the critical exponent 
for I is the same as that of the order parameter. 

This result can be understood in terms of the mean-field theory of fluctuations as follows. 
Consider the reduced pressure, n = (P — P c )/P c , close to the critical point. Expanding to 
first order in the reduced temperature, r = (T — T c )/T c , and to third order in the reduced 
density, r] — (p — p c )/p c , one finds: 



7T 



ar + brr] + cr] 3 + c'tt] 2 . (28) 
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Note that the conditions to find the critical point, 

On d 2 ir 



0, (29) 



drj drf 

have already been used to eliminate possible explicit linear and quadratic terms in 77. To find 
the pressure at phase coexistence, one can use Eq. (IT81) in Eq. fl28|) to find the symmetric 
gas and liquid coexistence points. One can then show that the coexistence pressure is given 
by: 

n\coex = aT. (30) 

Using the Clausius-Clapeyron relation, the latent heat close to the critical point becomes: 

n P 

l = ^f c (Pl-p 9 )(l+T)~(-Tf. (31) 

In this expression, we have considered that the difference in densities is the order parameter 
of the phase transition and that the latter is governed by the critical exponent, p\ — p g ~ 



(— t) 13 . As any mean- field theory |6|, |29|, self-consistent Hartree-Fock calculations of the 
nuclear matter liquid-gas phase transition yield a critical exponent (3 = 1/2. 

We have checked numerically that the critical exponent of I is indeed = 1/2. Fig. [7] 
shows a logarithmic plot of I versus the reduced critical temperature (symbols). Linear 
regression fits have been performed and are shown with lines. The correlation coefficients 
are close to 1 to within at least four digits in all cases. Such a linear behavior of the data 
confirms the scaling of / with r. The slopes of these lines have also been extracted and 
agreement with the derived from the coexistence line is good up to the third digit. To 
our knowledge, this is the first time that the critical exponent of the latent heat for nuclear 
matter is computed and that its equivalence to is confirmed numerically in the framework 
of a Hartree-Fock mean-field theory. However, one should keep in mind the limitations of 
the mean-field theory for the calculation of the critical exponents. 



V. SUMMARY AND CONCLUSIONS 



In this work, we have have analyzed in detail the latent heat of the liquid-gas phase 
transition of symmetric nuclear matter. The latent heat describes the amount of heat 
needed to transfer one nucleon from the liquid to the gas phase, and it can be used as 
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a further characterization of the phase transition. We have been motivated by experimental 
results which suggest that, in the phase transition occurring in nuclear multifragmentation 



collisions, a latent heat of around I ~ 2 — 8 MeV can be observed [4], [15|, H 

We have used self-consistent Hartree-Fock calculations at finite temperature to obtain 
numerical results. Four different characteristic mean-fields have been used to determine the 
inherent mean-field dependence of the results. The qualitative behavior of the latent heat as 
a function of temperature is very similar for all effective interactions. Within the Hartree- 
Fock approximation, the latent heats derived from different mean-field parametrizations fall 
within a narrow band when the temperature is scaled by T c and I is scaled by e$. 

In the T — > limit, the latent heat coincides with the binding energy per particle at 
saturation density. The latent heat can also be computed from the difference in entropies 
between the gas and the liquid phases. We have seen that the gas entropy dominates over the 
liquid one in a wide range of temperatures. At finite but low temperatures, the latent heat 
rises linearly with temperature. A careful analysis shows that the slope of this linear trend 
is 5/2, regardless of the mean-field parametrization or the many-body approximation. We 
have confirmed this trend with numerical finite temperature HF calculations. This model- 
independent result is valid for all normal fermionic systems. To our knowledge, it is the first 
time that this result is discussed. 

Knowing that a) I is positive and finite at T = 0, b) it has a positive slope near T = 
and c) it goes to zero at T = T c , necessarily implies the existence of a maximum in the 
latent heat. Mean-field numerical calculations suggest that this maximum is located around 
T ~ 0.6T C and that l max ~ 1.8eo- While the exact position of the well-defined maximum 
can depend on the effective force or on the many-body approach used to describe nuclear 
matter, its existence is guaranteed independently of the mean-field parametrization or the 
many-body approximation. Moreover, this result should be valid for any normal fermionic 
system that presents a liquid-gas phase transition. 

It is not particularly easy to understand physically why the latent heat should present a 
maximum as a function of temperature. One might interpret the presence of this maximum 
as a manifestation of the underlying NN interaction. At low temperatures, the latent heat 
rises because more work is needed to break the attractive part of the strong interaction 
when transferring a nucleon from the more structured liquid to the gas phase. At higher 
temperatures, the thermal motion of the particles would provide most of the work and thus 



21 



less latent heat is needed to break the bonding between particles in the liquid. Here, we 
have given an analytical demonstration of the existence of the peak. 

Concerning the behavior of the latent heat near the critical point, we have shown nu- 
merically that it can be described in terms of a single critical exponent. As predicted by 
the theory of phase transitions, this critical exponent is the same as the one associated to 
the order parameter. In the case of Hartree-Fock calculations in nuclear matter, (3 = 1/2, 
and there is a very good numerical agreement between both exponents. All in all, the la- 
tent heat can be characterized in the low, intermediate and (close to) critical temperature 
regimes from very basic principles. It is precisely this generic nature which might motivate 
the use of the latent heat as a tool in analyzing the liquid-gas phase transition in other 
normal fermionic systems. 

We are well aware that a connection between experimental observations of multifragmen- 
tation collisions and theoretical results of homogenous nuclear matter is not at all transpar- 
ent. Finite size effects play a capital role in determining the thermodynamical properties of 
nuclei and the latent heat is not an exception 36J. Nuclear matter values at its maximum 
suggest that the latent heat is up to 10 times higher than that extracted from different 



experimental analysis 
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19] . Even at the theoretical level, it is not clear how to define 

n 

a liquid and a gas phase in a self-confined system [9] . Nevertheless, one might naively expect 
that the appearance of a maximum in the latent heat might have a substantial suppression 
effect on the yields of light particles. We hope that the present study, where we have high- 
lighted model-independent and basic arguments for the latent heat, will encourage further 
experimental analysis of the latent heat in multifragmentation collisions. 
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